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Abstract. - An analytical result for the renormalization of the jump of the heat capacity AC/Cn by the 
anisotropy of the order p arameter is derived within the framework of the very recent model proposed by Posazhen- 
nikova, Dahm and Maki |cond-mat/0204272| submitted to Europhys. Lett.], for both oblate and prolate anisotropy. 



The graph of AC/Cn versus the ratio of the gaps on the equator and the pole of the Fermi surface, A e /A p , allows 
a direct determination of the gap anisotropy parameter A e /A p by fitting data from specific heat measurements 
AC/Cn- Using the experimental value AC/Cn = 0.82 ± 10% by Wang, Plackowski, and Junod [Physica C 355 
(2001) 179] we find A e /A p « 4.0. 



In a very recent e-print Posazhennikova, Dahm and Maki Q] discuss a model for the gap anisotropy in MgB2 , 
a material which has attracted a lot of attention from condensed matter physicists in the past two years. A 
central issue in this work JlJ is to propose an analytic model for analyzing thermodynamic behavior. Assuming 
a spherical Fermi surface, a simple gap anisotropy function is suggested, A(k) = A e /vT + Az 2 , where z = cos#, 
and 8 is the polar angle. This model leads to useful results for the temperature dependence of the upper critical 
held H C 2 and of the specific heat, which can be fitted to the experimental data, thereby determining the optimal 
anisotropy parameter A. Note that A — (A e /A p ) 2 — 1, with A p = A(z = 1) and A e = A(z = 0), and the gap 
ratio is parameterized as A e /A p = ^/l + A > 0. 

The aim of the present Comment is to provide a convenient analytical expression giving the possibility for 
determining A e /A p from the available data for the jump of the specific heat U. For the latter we derive the 
explicit formula 

Ep is the Fermi energy, k are the band energies, (1) is the density of states, £ is the Riemann zeta function, 
is analogous to Abrikosov's parameter (3a ||, and 12/7^(3) = 1.42613. . . is the sacramental BCS ratio. Then 
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Fig. 1 - Jump of the specific heat AC/Cn versus the "equatorial-to-polar" gap ratio A e /A p . For a given AC /Cn value 
we have oblate A e /A v > 1 and prolate A e /A p < 1 solutions. The solid line gives our present analytical result eq. (^) 
for the model by Posazhennikova, Dahm and Maki jjj]. The dash-dotted line is our analytical solution |H) for the model 
by Haas and Maki [^|. The dashed line is the jump ratio AC/Cn = 0.82 ± 10% measured by Wang, Plackowski, and 
Junod p|, with the shaded area showing the experimental error bar. 




following the weak-coupling BCS approach [|l],y we derived the explicit analytic expressions valid for A > 0, 
and — 1 < A < 0, respectively 



AC (A) 



12 2(l + fr 2 )(arctanb) 
7C(3) b 2 + &(l + & 2 )arctan&' 

12 2(1 -p 2 ) (tanbT 1 ?}) 2 
7C(3)p 2 +p(l-p 2 )tanh- i p' 




(2) 



For a given specific heat jump, this expression leads to two solutions (oblate, A e /A p > 1, and prolate, A e /A p < 
1). The relevant example is shown in fig. [l]; the function AC(A)/Cn is tabulated in ref. Q. The analysis of the 
angular dependence of H c i [^[0] performed in the commented paper Q| unambiguously demonstrates that one 
has to analyze only the "oblate" case. Thereby the experimentally reported value in ref. Q AC/Cn — 0.82±10% 
gives A ss 16 and A e /A p «6« 4.0 ± 10%. For this significant anisotropy, the "distribution" of Cooper pairs 
(|A k | ) <* !/[*£+ (kF/b) 2 } has a maximum at k z = 0. This general qualitative conclusion is in agreement with 
the hints from band calculations that the maximal order parameter is concentrated in an almost two-dimensional 
electron band, but all bands Eb,k take part in the normal specific heat per unit cell Cn = i^ 2 /3)kgT(l) . 

For the two-band model, advocated for the first time for MgB2 in ref. M, eq. (]]]) gives (to within a typo- 
graphical correction) the result by Moskalenko |J 

C N 7C(3) (pi+p 2 )(|A 1 |Vi + |A 2 | 4 P2) z±x+{l~x)' A 2 pi + p 2 W 

and pi and p 2 are the densities of states for the two bands. Taking for an illustration x — 0.515 and AC/Cn = 
0.82, eq. @ gives Ai/A 2 ~ 4.0 in agreement with A e /A p « 4.0 obtained using eq. (||). Thus the gap ratios are 
model-independent. For a survey on a set of parameters see Table I in ref. jlO). Certainly the jump of the heat 
capacity alone cannot be an arbiter for the validity of any model, so subtleties, e.g., related to strong coupling 
effects and other anisotropies, can be hidden in the parameters spread in the table mentioned. 
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